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Leader~mRNA Junction Sequences Are Unique for Each Subgenomic mRNA Species 


in the Bovine Coronavirus and Remain So Throughout Persistent Infection’ 
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The common leader sequence on bovine coronavirus subgenomic mRNAs and genome was determined. To examine 
leader-mRNA junction sequences on subgenomic mRNAs, specific oligadeoxynucleotide sets were used in a polymer- 
ase chain reaction to amplify junction sequences from either the positive-strand MRNA (eight of nine total identified 
species) or the negative-strand anti-mRNA (six of the nine species), and sequenced. The mRNA species studied were 
those for the N, M, S, and HE structural proteins and the 9.5-, 12.7-, 4.8-, and 4.9-kDa putative nonstructural proteins. 
By defining the leader-mRNA junction sequence as the sequence between (i} the point of mismatch between the leader 
and genome and (ji) the 3’ end of the consensus heptameric intergenic sequence [(U/A)C(U/C)AAAC)], or its variant, a 
unique junction sequence was found for each subgenomic mRNA species studied. In one instance (mRNA for the 
12.7-kDa protein) the predicted intergenic sequence UCCAAAC was not part of the junction region, and in its place was 
the nonconforming sequence GGTAGAC that occurs just 15 nt downstream in the genome. Leader-mRNA junction 
sequences found after 296 days of persistent infection were the same as those found during acute infection (<18 hr 
postinfection). These data indicate that, in contrast to the closely related mouse hepatitis virus, the bovine coronavirus 
maintains a stable leader-mRNA junction sequence for each mRNA. Interestingly, this stability may be related to the 
fact that a UCUAA sequence element, postulated by others to be a regulator of the leader-mRNA fusion event, occurs 
only once within the 3’ flanking sequence of the genomic leader donor and once at intergenic sites in the bovine 


coronavirus genome, whereas it occurs two to four times at these sites in the mouse hepatitis coronavirus. 
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INTRODUCTION 


It is well established that the 3’ coterminal subgeno- 
mic mRNAs of coronaviruses have a common 5’ termi- 
nal leader sequence of approximately 65 to 90 nucleo- 
tides (the length depends on the species of corona- 
virus) which is not encoded by a colinear region on the 
genome, but rather by the 5’ end of the genome (for 
reviews see Lai, 1990; and Spaan et a/., 1988). The 
leader on each subgenomic mRNA makes up only a 
portion of the 5’ untranslated region (Fig. 1) which, in 
turn, differs in length for each mRNA species (Table 1). 
Assuming that subgenomic mRNA molecules are de- 
sived from full-length genome during virus replication, 
then a mechanism must exist whereby a copy of the 
genomic leader becomes fused within the 5’ untrans- 
lated region to the subgenomic MRNAs during virus 
replication. For each of severat models that have been 
proposed to explain the fusion event (described be- 
low), a conserved heptameric intergenic sequence oc- 


' Sequence data from this articie have been deposited with Gen- 
Bank Data Library under Accession Number M62375. 
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curring in the genome at varying distances upstream 
from the start cadon of the various genes is thought to 
be an important sequence element determining the fu- 
sion site, since the point of fusion in the 5’ untranslated 
region of the mRNA is within or just upstream of the 
consensus intergenic sequence. 

The three major hypotheses that have been put for- 
ward to explain the leader-mRNA fusion process are 
the following: () The discontinuous transcription hy- 
pothesis states that the polymerase discontinues syn- 
thesis after transcribing leader from the 3' end of the 
minus strand, and jumps to a new transcription site at 
conserved intergenic heptameric sequences on the 
antigenome to continue synthesis of the mRNA (Baric 
et a/., 1983; Lai, 1988; Spaan et a/., 1983). This model 
has been slightly modified to become the leader-prim- 
ing hypothesis which states that teader is made in ex- 
cess, becomes free, and serves to prime transcription 
at the intergenic sites on the minus strand in cis or in 
trans (Makino et a/., 1986). (ii) The minus-strand splic- 
ing hypothesis states that full-length minus-strand 
anti-genomes become spliced either in cis or in trans to 
mRNA-length minus-strand molecules that then serve 
as templates for synthesis of mRNAs (Sawicki and Sa- 
wicki, 1990). (ii) The plus-strand splicing hypothesis 
states that progeny genomes become directly spliced 
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Fig. 1. Schematic definition of the leader-mRNA junction se- 
quence. (A) Postulated ‘‘free’’ leader containing the common leader 
sequence (65 nt for BCV, identified in bold type), the flanking con- 
sensus intergenic sequence, and extended sequence of unknown 
tength. (B) Alignment of leader and remaining mRNA with the anti- 
genome after synthesis of mRNA by the postulated leader-primed 
mechanism of transcription. BCV N mRNA is used here for illustra- 
tion. In this case, the teader-mRNA junction sequence is AU- 
CUAAAC. (C) 5’ end of the completed BCV N mRNA showing the 
leader sequence and the leader-mRNA junction sequence as com- 
ponents of the 5’ untranslated region. 


either in cis or in trans into subgenomic MRNA mole- 
cules (Baric et a/., 1983; Lai, 1988). Of these hypothe- 
ses, the leader-priming hypothesis has been used the 
most widely to explain the origin of leader-containing 
subgenomic mRNA molecules. With any one of the 
three models, however, it might be expected that the 
fusion event would be an imprecise one and that vari- 
ant fusion junctions for any given mRNA species would 
be observed from a collection of imprecise fusion 
events. For the leader priming model, for example, it 
could be postulated that the free leader might misalign 
during the priming process and create a different 
leader-mRNA junction. For the splicing models, these 
variations could arise from alternative splice sites. 
Recent studies on the mouse hepatitis virus (MHV} 
have shown that multiple leader-mRNA junction se- 
quences do in fact exist for certain of the subgenomic 
mRNA species (Makino et a/., 1988). These findings 
were interpreted in light of the leader-priming model 
and explained by the existence of muitiple 5°UCUAA 
sequence sets in the putative free leader that align in 
optional ways with multiple 3'AGAUU repeats occur- 
ring within intergenic regions on the antigenome (Ma- 
kino and Lai, 1989; Baker and Lai, 1990). Thus, the 
resulting mRNA species had two, three, or four UCUAA 
repeats within the leader-mRNA junction sequence. It 
was further documented for the JHM strain of MHV that 
transcription of one subgenomic mRNA species, 


mRNA 2b (the transcript for the HE protein), occurred 
when two UCUAA repeats existed in the genomic 
leader (and thus also the presumed free leader), but 
not when three existed (Shieh et a/., 1989). Other varia- 
tions of this phenomenon have been described for 
MHV (LaMonica et a/., 1992). The number of UCUAA 
repeats was thus postulated to be part of amechanism 
for regulating the rate of transcription initiation (La- 
Monica et a/., 1992; Shieh et a/., 1989). 

We have sought to study the leader-mRNA junc- 
tions on bovine coronavirus (BCV) mRNAs for several 
reasons. (i) Sequence analyses may provide clues to 
the leader fusion mechanism. (ii) BCV is a close relative 
of MHV, on the basis of amino acid sequence similari- 
ties in the proteins (Abraham et a/., 1990a,b; Kienzle et 
al., 1990; Lapps et a/., 1987), and a comparative analy- 
sis of leader and leader-mRNA junction sequences 
may point out important conserved structural features 
that might suggest leader function(s). This is especially 
important with regard to fusion sequence variations 
and the regulation of gene expression. We notice, for 
example, no variation in the levels of HE protein or HE 
MRNA (species 2b) production in different subclones 
or passages of BCV (Hogue et a/., 1989; Hofmann et 
al., 1990; data not shown). During BCV replication both 
are invariably expressed at high levels. (iii) Highly vari- 
able junction sequences would challenge the notion of 
MRNA replication (Sethna et a/., 1989). We and others 
have recently shown that coronavirus subgenomic 
mRNAs have minus-strand counterparts that appear to 
be active in transcriptional or replicational complexes 


TABLE 1 


PROPERTY OF 5'-UNTRANSLATED REGIONS AND LEADER-mRNA JuNc- 
TION SEQUENCES ON BOVINE CORONAVIRUS GENOME AND SUBGENOMIC 
mRNAs? 


mRNA Length of Leader-mRNA 
species S'UTR (nt) junction sequence 
Genome 210 -TAAAC- 
32K ND ND 
HE 79 -CTAAAC- 
$ 70 -ATAATCTAAAC- 
(4.9K)? (no start codon) (-CTAAGT-} 
4.8K 75 -TAAAC- 
12.7K 124 -TAGAC- 
9.5K 193 -AATCCAAAC- 
M 73 -TAATCCAAAC- 
N 77 -ATCTAAAC- 


*The 5’ leader sequence of 65 nt is identical on all sequenced 
BCV RNA species. The leader-mRNA junction sequence is defined 
as that sequence which lies between (j) the point of mismatch be- 
tween leader and genome, and (ii) the 3’ end of the consensus inter- 
genic sequence or its variant. 

» Although the mRNA species for the 4.9-kDa protein was found, it 
is apparently nonfunctional since it contains no start codon (as ex- 
plained in the text). ND, not determined. 
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MHV-A59 «5 -uaUaaGAGUGRUU-gGCGUcCGUal gUalCCuCUcaACucUagAaCUCUUG-UAGUUuaaa-UcURA-----| WCURARCUUUARGGAUG, .. 
MHV-JHM 5‘ -uaUaaGAGuGAUU-g6CGUcCGUal gUalCCuCUcuACucUaaAaCUCUUG-UAGUUUada-UcUAAUcuaaUCUARRCUUUAAGGAUG. .. 
HCV-0C43 5 --UUGUGRGCGRaguUGCGUGCGUGCAU-CCCGECUUCAC-cU-GAUCUCUUGUUAGAUCUUUUR cUAAMcuaa)/CURARUUUURAGGAUG. .. 
BCV 5’ GAUUGUGAGCGRUU-UGCGUGCGUGCAU-CCCGCUUCAC--U-GAUCUCUUGUURGRUCUUUUUAURA-----| CUAARCUUUARGGAUG. . . 
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Fig. 2. Alignment of the leader-containing 5’ untranslated regions on the N mRNAs for the bovine coronavirus and three of its close relatives. 
Numbers on top refer to the composite sequence, and on the bottom to the BCV sequence. Nucleotides that are identical to the BCV sequences 
are shown in upper case letters. The putatively critical UCUAA sequences are boxed, and the N start codons are underlined. (One apparently 
noncritical UCUAA set beginning at base position 39 in the MHV strains is not boxed.) The first 65 nt of the BCV sequence make up the BCV 
leader. The first 74 nt of the illustrated BCV sequence are identical to the first 74 nt on the BCV genome. 


in the infected cell (Hofmann et a/., 1990; Sawicki and 
Sawicki, 1990; Sethna et a/, 1989). If subgenomic 
mRNAs do undergo replication, then this would be one 
mechanism by which MRNAs could maintain a stable 
structure during persistent infection. Conversely, varia- 
tions in leader-mRNA junction sequences over tire 
would be difficult to explain if mRNAs were perpetu- 
ated only by a replication mechanism. 

In experiments described here, the common leader 
sequence on the BCV subgenomic mRNAs and ge- 
nome was determined, and by sequencing PCR ampli- 
fied junctions, leader-mRNA junction sequences for 
eight of nine species of BCV subgenomic mRNAs were 
determined. The junction sequences on the anti- 
mRNAs (six of nine species) were also determined. We 
learned that the leader-mRNA junction sequence was 
unique for each subgenomic mRNA and remained so 
throughout persistent infection. The intergenic can- 
sensus sequence (U/A)C(U/C)AAAC occurs no more 
than once in any given BCV leader—mRNA junction se- 
quence, and on one mRNA species does not occur at 
all. The leader-mRNA junction sequences were there- 
fore remarkably stable in the bovine coronavirus and, 
by comparison with published data on the mouse coro- 
navirus, this stability appears to be related to the singu- 
lar occurrence of the UCUAA sequence element in the 
genomic leader 3’ flanking sequence and at each inter- 
genic region in the genome. 


MATERIALS AND METHODS 


Virus and cells and preparation of RNA from 
infected cells 


The Mebus strain of BCV was grown on human rec- 
tal tumor (HRT) cells as previously described (King and 
Brian, 1982; Lapps et a/., 1987). A persistent infection 
in cell culture was established, and RNA was extracted 
from acutely and persistently infected cells as previ- 
ously described (Hofmann ef a/., 1990). In the current 
study, persistently infected cells had been carried for 
74 passages (296 days). 


Determination of the BCV leader sequence on the 
N, M, and S mRNAs and genomic RNA 


To establish the BCV leader sequence on N mRNA 
(Fig. 2), oligodeoxynucleotide 1, which binds to a re- 


gion beginning 34 nt downstream from the N transla- 
tion initiation codon (Lapps et a/., 1987) (Fig. 3), was 5’ 
end-labeled with polynucleotide kinase and [y-@PJATP 
(>3000 Ci/mmol; ICN Pharmaceuticals), and used as 
primer for an extension reaction on cytoplasmic RNA 
harvested at 6 hr postinfection (Hofmann et a/., 1990). 
The extended primer was isolated by electrophoresis 
on a 12% polyacrylamide sequencing gel and se- 
quenced by the chemical method (Maxam and Gilbert, 
1980}. This is essentially the same method used to 
determine the leader sequence of MHV (Lai et a/., 
1984), the avian infectious bronchitis virus (IBV) (Brown 
et a/., 1984), and the porcine transmissible gastroen- 
teritis virus (TGEV) (Sethna et a/., 1991). To confirm the 
leader sequence, the experiment was done with ex- 
tended products from the M mRNA template (Lapps et 
al., 1987). For this, oligodeoxynucleotide 2 which 
binds to a region beginning 19 nt downstream from the 
predicted M start codon (Fig. 3) was likewise used. 

Because it was impossible with this experimental 
approach to unequivocally resolve the 5’ terminal 
bases, a method was developed in which extended 
primers in separate experiments for each the N, M, and 
S mRNAs were ligated into head-to-tail multimers 
(dimers and trimers) with RNA ligase, and PCR-ampli- 
fied molecules representing the joined regions were 
cloned and sequenced by the dideoxynucleotidyl 
chain-termination method (Hofmann and Brian, 
1991a,b). 

To determine the leader sequence on the BCV ge- 
nome, the del site-containing oligonucleotide 5’CCT- 
CCAAATCATATGGACGTGTATTC3’ which binds to a 
region on genomic RNA beginning 456 nt downstream 
from the 9’ end, and the EcoRI, Kpnl site-containing 
oligonucleotide5’'CGGAATT CGGTACCGATTGTGAGC- 
GATTTGCGTGCG3' which binds to the first 22 nt of the 
antileader, were used in a PCR to generate a 495 nt 
fragment that was cloned and sequenced. Genomic 
RNA was obtained from a virus purified from the super- 
natant fluids of persistently-infected cells at 120 days 
postinfection (Hofmann et a/., 1990). At this passage 
level, there was no detectable defective interfering 
RNA species in the cells or virus. The 5’ end of the 
genome was originally obtained from a cloned and se- 
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quenced defective-interfering RNA which showed a 5’ 
terminal nucleotide sequence identity of 64% with the 
5’ terminal 498 nt of the MHV genome (Chang et @/., 
submitted for publication; Soe et a/., 1987) and it was 
from this sequence that the Ndel site-containing oligo- 
nucleotide (described above) was obtained. The 5’ ter- 
minal 455 sequence of the genome is identical to that 
at the defective-interfering RNA. Since some of the 5’ 
23 nt of the genomic leader may have been dictated by 
the antileader primer, this portion of the leader remains 
to be confirmed. 


Asymmetric PCR amplification and sequencing of 
leader-mRNA junction sequences 


For analysis of leader-mRNA junction sequences, 
oligodeoxynucleotides 1 and 2 (described above), and 
oligodeoxynucleotides 3 through 8 (described in Fig. 3) 
which bind respectively to regions within the open 
reading frames of the 9.5-, 12.7-, 4.8-, 4.9-kDa, S, and 
HE proteins, (Abraham et a/., 1990a,b; Kienzle et a/., 
1990) were prepared and used for first-strand cDNA 
synthesis and for subsequent thermocycling reac- 
tions. A 26-mer oligodeoxynucleotide (5’°GAGCGATT T- 
GCGTGCGTGCATCCCGC3’) which binds to nt 7 
through 32 of the antileader (counting from the 3’ end 
of the antileader) was prepared and used as the sec- 
ond primer in all thermocycling reactions. 

For first-strand cDNA synthesis from mRNA tem- 
plate, 75 wg of total cytoplasmic RNA (stored as an 
ethanol precipitate) was pelleted, washed with 80% 
ethanol, dried, dissolved in 12.25 ul water, denatured 
by heating at 70° for 5 min, and quick-cooled on ice. 
First-strand cDNA was synthesized in a 30-p! reaction 
mixture containing 12.25 ul (75 ug) RNA, 5 ul (~50 
pmol) first-strand mRNA-specific primer, 3 ul 10x re- 
verse transcriptase buffer (1X reverse transcriptase 
buffer = 50 mM Tris-HCl, pH 8.3, at 42°, 7 MM MgCl, 
40 mM KCl, 1 mM DTT, and 0.1 mg BSA per ml), 6 ul of 
the four dNTPs {at 5 mM each), 0.75 ul (30 units) RNa- 
sin (Promega), and 3 ul (24 units) AMV reverse tran- 
scriptase (Promega), for 2 hr at 42°. The reverse tran- 
scriptase reaction was stopped by incubation at 75° 
for 15 min, and the RNA was removed by adding 2 yl 
RNase A (at 2 mg per ml; Sigma) and incubating the 
mixture at 37° for 1 hr. ssDNA was extracted once with 
phenol chloroform and once with chloroform. 

For second-strand cDNA synthesis and PCR ampilifi- 
cation, a 50-ul PCR reaction mix containing 25.75 ul 
H,0, 5 wl 10X PCR buffer (1X PCR buffer = 10 mM 
Tris-HCl, pH 8.3, 50 mM KCI, 1.5 mM MgCl, 0.01% 
gelatin, and 0.01% Triton X-100), 8 wl dNTPs (at 1.25 
mM each), 5 ul of the antileader-specific primer (at 10 
pM), 5 wl of the respective mRNA-specific primer (at 10 
uM), 1 ul of first-strand synthesis mix, and 0.25 ul (1.25 
units) Tag DNA polymerase (Promega) was prepared 


and incubated for 30 cycles in the thermal cycler, 1 min 
at 94°, 1 min at 55°, and 1 min at 72° for each cycle. 

For analysis of minus-strand RNA, the same proce- 
dures as described above for first and second-strand 
cDNA synthesis and PCR amplification were used, ex- 
cept the order of primer usage was reversed. 

To examine the sizes of the amplified doubte- 
stranded DNA products and to obtain a small amount 
of purified product for asymmetric PCR amplification, 
20 yl of the PCR reaction mix was electrophoresed ina 
4% agarose gel (NuSieve, FMC) containing Tris—borate 
buffer and visualized by ethidium bromide staining. The 
predicted sizes of the amplified products were 129 bp 
for the N mRNA, 104 bp for M mRNA, 235 bp for the 
9.6-kDa protein mRNA, 140 bp for the 12.7-kDa pro- 
tein MRNA, 96 bp for the 4.9-kDa protein mRNA, 99 bp 
for the 4.9-kDa protein mRNA, 106 bp for the S mRNA, 
and 124 bp for the HE mRNA. 

For asymmetric amplification of DNA by PCR in prep- 
aration for DNA sequencing, a modified method of 
Innis et a/., (1988) was used. To a 50-ul reaction mix 
containing 42.2 yl water, 5 ul 10X PCR buffer, 0.8 ul 
dNTP's (at 1.25 mM each), 1 ul of the antileader-speci- 
fic primer (the limiting primer; 0.1 #4), and 1 wl of the 
mRNA-specific primer (10 u//) was added a small inoc- 
ulum of amplified DNA from the agarose gel. The DNA 
was transferred by stabbing the band with a micropipet 
tip and rinsing the tip in the asymmetric reaction mix. 
The reaction mix was overlaid with mineral oi and in- 
cubated for 40 cycles in a thermocycler, 30 sec at 94°, 
1 min at 55°, and t min at 72°. The asymmetrically 
amplified single-stranded DNA was separated from un- 
incorporated dNTPs by chromatography in a water- 
equilibrated Sephadex G25 Spin column, dried, and 
redissolved in 20 wl 1x PCR buffer. 

For sequencing asymmetrically amplified DNA, the 
protocol of Innis et a/. (1988; Hofmann and Brian, 
1991b) for dideoxynucleotidy!] DNA sequencing with 
Taq DNA polymerase was used. The 5’ end of the 
leader-specific primer was labeled with °2P using poly- 
nucleotide kinase and [y*2PJATP (+3000 Ci/mmol: 
ICN) and used in the extension reactions. DNA was 
electrophoresed on 8% polyacrylamide gels contain- 
ing 50% urea. 


RESULTS 
BCV leader sequence on mRNAs and genome 


The entire 77-nt untranslated region on the BCV N 
mRNA is shown in Fig. 2. A majority of the BCV 5’ 
untranslated region for the N and the M mRNAs was 
determined by chemical sequencing, but the 5’ termi- 
nal five nucleotides for these mRNA species as well as 
for the S mRNA were established by sequencing li- 
gated head-to-tail cDNA products prepared from 
mRNA obtained during acute infection (Hofmann and 
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Brian, 1991a,b). From these data, and from the leader— 
mRNA sequence data on all but one of the subgenomic 
mRNAs described below, we conclude that the BCV 
common leader on subgenomic mRNAs is the first 65 
nucleotides of the BCV sequence shown in Fig. 2. Fur- 
thermore, by sequencing a cloned 495-nt PCR frag- 
ment of the 5’ end of the genome, it was established 
that the genomic 5’ terminus is identical to the N mRNA 
5’ terminus for nt 23 through 74, and probably also 
identical for nt 1-22 (the uncertainty arises because 
this is the region to which the primer bound) (Fig. 2). 

Since the leader-mRNA fusion sites appear to be 
within or just upstream of a singular consensus hepta- 
meric intergenic sequence [(U/A)C(U/C)AAAC)], the 
leader-mRNA junction sequence was defined in this 
study as that region extending from the point of diver- 
gence between the leader and genome, to the 3’ end of 
the heptameric intergenic sequence (or its variant) 
(Fig. 1). 


Leader-mRNA junction sequences are unique for 
each subgenomic mRNA species 


The results of each sequenced junction for eight of 
nine species of subgenomic MRNA are shown in Fig. 3 
and Table 1. The junction sequence on the subgeno- 
mic mRNA of the 32-kDa nonstructural protein de- 
scribed by Cox et a/. (1989) was not studied since the 
sequence of this gene in the Mebus strain of BCV had 
not been determined. 

From our analyses, the following points emerge: (i) 
The leader—junction regions range in length from 5 nt 
for the MRNAs of the 12.7- and 4.8-kDa proteins to 11 
nt for the S mRNA, and the leader—junction sequence 
is unique for each of the eight subgenomic mRNAs 
studied. (ii) The positions of five of the leader-mRNA 
junction Sequences were as predicted from the posi- 
tions of a (U/A)C(U/C)AAAC consensus sequence on 
the genome, from the proximity of these sequences to 
the putative AUG initiator codons on the deduced 
genes, and from the size of the mRNAs as judged by 
Northern analyses. These were for the mRNAs of the 
N, M, 9.5-kDa, S, and HE proteins (Abraham et a/., 
1990a,b; Kienzle et a/., 1990; Lapps et a/., 1987). (iii) 
Three of the eight junction regions, however, were dif- 
ferent from those predicted. (a) The mRNA for the 12.7- 
kDa protein was found to be joined to the leader se- 
quence 21 nt downstream from the predicted 
UCCAAAC sequence (at boxed position 5953 in Fig. 3), 
and joined at the unusual consensus variant GGUA- 
GAC (Fig. 4). The resulting mRNA would be expected 
to be functional since the first AUG on the sequence 
predicts the start of the identified 12.7-kDa protein. (b) 
The mRNA having the 4.8-kDa ORF as the upstream 
ORF was not predicted earlier since no (U/A)C(U/ 
C)AAAC consensus sequence was found near the pu- 
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Fig. 4. The leader-mRNA junction sequence for the mRNA of the 
12.7-kDa protein. The asymmetrically amplified DNA was se- 
quenced and the leader-mRNA junction sequence is shown. 


tative AUG initiation codon, and the 4.8- and 4.9-kDa 
ORFs together appeared to be the result of a point 
mutation disrupting an ancestral 11-kDa ORF of the 
kind found in MHV (Abraham et a/., 1990). The 4.8-kDa 
ORF junction region appears to utilize AGUAAAC as 
the intergenic sequence which requires a mismatch in 
the second nucleotide of the consensus sequence. (c} 
Most surprising is an intergenic sequence near the be- 
ginning of the AUG for the 4.9-kDa protein MRNA. This 
junction did not utilize the predicted ACCAAAC con- 
sensus sequence beginning 323 nt upstream from the 
AUG, but rather used the ACUAAGU sequence begin- 
ning 4 nt downstream from the putative AUG. This 
transcript is therefore probably not used for protein syn- 
thesis since the first AUG in the transcript, even though 
itis ina good sequence context for initiation of transla- 
tion (Kozak, 1991), is 102 nt downstream and would 
yield a peptide of only four amino acids. 
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To confirm the existence of the leader—-mRNA junc- 
tion sequences, oligodeoxynucleotides that would 
prime and amplify the minus-strand copy of each 
leader-mRNA junction (Hofmann et a/., 1990) were 
used and the asymmetrically amplified DNA was se- 
quenced. Abundant discrete PCR products were ob- 
tained from the minus strands for all but the 4.9- and 
4.8-kDa protein mRNAs, and hence no DNA sequence 
was obtained from these two species. These results 
are consistent with our earlier RNA blotting analysis 
(Hofmann et a/., 1990), and suggest that minus-strand 
copies of the specific 4.9- and 4.8-kDa mRNA species 
either do not exist or are present in very small numbers. 


Leader-mRNA junction sequences remain 
unchanged after 296 days of persistent infection 


To examine the leader-mRNA junction sequences 
after a long period of persistent infection, cytoplasmic 
RNA extracted from cells at 296 days (74 passages) 
postinfection was used in an identical manner to that 
described above to separately identify junction se- 
quences on plus- and minus-strand subgenomic RNA 
species. Except for our failure to confirm sequences on 
minus-strand copies of mRNAs for the 4.9- and 4.8- 
kDa proteins, junction sequences were found to be 
identical to those described in Fig. 3. 


DISCUSSION 


We have established the sequence of the 5’ untrans- 
lated region for the bovine coronavirus N and M sub- 
genomic mRNAs by using a combination of three sepa- 
rate methods, direct chemical sequencing of an ex- 
tended primer, sequencing of a cloned extended 
primer, and dideoxynucleotidyl sequencing of DNA 
asymmetrically amplified from an extended primer. By 
using the latter two methods, the 5’ untranslated re- 
gion on the S mRNA was alsa determined (Hofmann et 
al, submitted for publication; and Fig. 3). From these 
data, we conclude that the common leader on BCV 
genome and subgenomic mRNAs is 65 nt in length 
(Fig. 2). The BCV leader is now the seventh coronavirus 
leader to be sequenced, and its sequence shows more 
similarity to the leaders of the antigenically closely re- 
lated human coronavirus HCV-OC43 (Kamahora et a/., 
1989), MHV-JHM (Shieh et a/., 1987), and MHV-A59 
(Lai et af, 1984), than to those of the more distantly 
related porcine transmissible gastroenteritis virus 
(Sethna et a/., 1991), avian infectious bronchitis virus 
(Brown et a/., 1984) and HCV-229E (Schreiber et a/., 
1989). The leader and flanking consensus sequences 
of the BCV-related coronavirus species can be com- 
pared by aligning the 5’ untranslated regions of their 
respective N mRNAs (Fig. 2). From this it can be seen 
that the BCV leader has more identity with HCV-OC43 
(85%) than with MHV-A59 or MHV-JHM (67 and 63%, 


respectively), suggesting that BCV and HCV-OC43 
share a common branch point during the evolution of 
this antigenic virus cluster. One notable feature is that 
whereas the pentanucieotide UCUAA is found as a se- 
quentially repeated sequence set near the 3’ end of the 
leaders of HCV-OC43 (three times), MHV-JHM (three 
times), and MHV-A53 (two times), it is present only 
once near the leader 3’ end on BCV mRNAs and ge- 
nome (the UAA portion actually flanks the leader) 
(Fig. 2). 

Our analyses of the bovine coronavirus feader— 
mRNA junction sequences reveal first of all that a 
unique junction sequence exists for each MRNA, and 
secondly that the uniqueness is maintained through- 
out at least 296 days of persistent infection. Regard- 
less of whether coronavirus mRNAs arise by only con- 
tinuous de novo synthesis (transcription) from a Jarger- 
than-unit size template (Lai, 1988), or by a two stage 
process in which the leader-mRNA fusion event is fol- 
lowed by faithful mRNA replication (Hofmann et a/., 
1990; Sethna ef a/., 1989), the remarkable stability (.e., 
invariance) of the BCV leader—-mRNA junction se- 
quences, in our view, is probably mechanistically re- 
lated to the singular occurrence of the UCUAA se- 
quence element in the putative leader donor sequence 
and at intergenic sites in the genome. This conclusion 
is based on (i} a comparison of BCV mRNA structure 
with the variable nature of leader-mRNA junction se- 
quences in MHV which correlates with multiple 
UCUAA sequence sets in the putative free leader (do- 
nor) and at intergenic sites (Baker and Lai, 1990; Ma- 
kino et a/., 1988; Makino and Lai, 1989; Monica et a/., 
1992; Shieh et a/., 1989), and {ii) recent in vitro muta- 
genesis studies in which an altered UCUAA intergenic 
sequence disrupted formation of subgenomic tran- 
scripts from a replicating defective-interfering RNA 
template Joo and Makino, 1992; Makino et a/., 1991). 
While the stability of the BCV leader-mRNA junction 
sequences in and of itself is consistent with the notion 
of mRNA replication, two other observations leave 
open the possibility that mRNAs arise only through a 
continuous transcriptional process requiring an accom- 
panying leader fusion step. {i) The stability of the BCV 
leader—-mRNA junction sequences contrasts sharply 
with a pattern of hypervariability of the 5’ end of sub- 
genomic mRNAs observed throughout establishment 
and maintenance of persistent infection (Hofmann et 
al., submitted for publication}. One would expect that 
replication should faithfully copy the termini as well as 
internal (leader-mRNA) junction sequences. (ii) Direct 
evidence to date does not support the notion that BCV 
mRNA molecules undergo complete replication. 
Whereas a cloned subgenomic defective RNA mole- 
cule carrying a reporter sequence undergoes replica- 
tion after transfection into infected cells, copies of 
cloned mRNA carrying the same reporter sequence do 
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not (Chang ef a/., submitted for publication). Since 
minus-strand anti-mRNA molecules exist, they may 
serve as templates for transcription by the use of an 
internal promoter, however (in the same manner as an- 
tigenome), but in this case strict control of the fusion 
event would likewise be expected. 

Although we have not clarified the mechanism of 
leader fusion onto mRNAs in these studies, two fea- 
tures of our data are difficult to reconcile with a leader- 
priming mechanism that depends on base pairing be- 
tween the complementary heptameric consensus se- 
quences on plus (free leader) and minus (antigenome 
or antimRNA) strands. The first is that the most favor- 
able heptameric intergenic sequence upstream of the 
12.7-kDa protein gene, UCCAAAC (which would allow 
a pairing of six bases), appears to have been bypassed 
in favor of the unusual GGUAGAC intergenic sequence 
(which allows a pairing of only four bases), occurring 
just 15 nucleotides downstream, for becoming part of 
the leader-mRNA junction. More than base pairing ap- 
pears to have directed this alignment. The second fea- 
ure difficult to reconcile is that the leader-mRNA junc- 
tion sequences we find are not dictated by the putative 
free leader but rather by the genomic sequence. The 
leader priming hypothesis postulates that the inter- 
genic consensus sequence on mRNAs will be derived 
from the free leader primer since the primer is ex- 
ended during the transcription process. From our stud- 
ies we learn that rather than the expected TAAAC be- 
coming part of the leader-mRNA junction sequence 
as predicted from the putative free leader sequence), 
CAAAC becomes part of the junction on the M and 
9.5-kDa mRNAs, TAGAC on the 12.7-kDa mRNA, and 
TAAGT on the 4.9-kDa mRNA. A similar problem has 
also been noted for MHV (Makino et a/., 1988) which 
led Makino et a/. to postulate a corrective proof-reading 
enzymatic action as part of the coronavirus transcrip- 
tase. In addition, a corollary of the base-pairing-based 
leader-priming hypothesis that states that the rate of 
mRNA synthesis will be directly proportional to the de- 
gree of base pairing at the intergenic consensus region 
(Budzilowicz et a/., 1985; Shieh et a/, 1987), is not 
fulfilled for BCV. Whereas the abundance of BCV 
mRNAs at the peak time of MRNA synthesis was mea- 
sured to be inversely proportional to mRNA length 
{Hofmann et a/., 1990), namely N > M > 9.5-kDa > 
12.7-kDa > S > HE, the degree of base pairing found 
within the leader-mRNA junction regions in this study 
predicts the relative mRNA abundances to be S > M> 
9.5-kDa > N > 4.9-kDa = HE > 12.7-kDa = 4.8-kDa 
(Fig. 3). Base pairing outside the region of the leader- 
mRNA junction, however, may contribute to the puta- 
tive priming event and these have not been considered 
in this prediction. 

Although, in our view, the mechanism of leader fu- 
sion onto MRNA remains to be established, our data 


lend support to the idea that the UCUAA sequence is 
important to the fusion process, and that the fusion 
event is strictly controlled when only a singular UCUAA 
sequence element is present near the fusion site either 
in the 3’ flanking sequence of the putative leader donor 
or at the intergenic template. The bovine coronavirus 
appears to be a useful system with which to examine 
the details of this controlled event since the UCUAA 
set is present only once at these sites in this corona- 
virus. 
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